Applications of nuclear magnetic resonance spectroscopy: from drug discovery to protein structure and dynamics.

The versatility of nuclear magnetic resonance (NMR) spectroscopy is apparent when presented with diverse applications to which it can contribute. Here, NMR is used i) as a screening/ validation tool for a drug discovery program targeting the Phosphatase of Regenerating Liver 3 (PRL3), ii) to characterize the conformational heterogeneity of p53 regulator, Murine Double Minute X (MDMX), and iii) to characterize the solution dynamics of guanosine monophosphate kinase (GMPK). Mounting evidence suggesting roles for PRL3 in oncogenesis and metastasis has catapulted it into prominence as a cancer drug target. Yet, despite significant efforts, there are no PRL3 small molecule inhibitors currently in clinical trials. This work combines screening of an FDA-approved drug panel and the identification of binders by protein-observed NMR. FDA-approved drugs salirasib and candesartan were identified as potent inhibitors in in vitro inhibition and migration assays while a weak inhibitor, olsalazine, was identified by NMR as the first small molecule inhibitor to directly bind PRL3. NMR was also used to validate the binding of additional compounds identified as experimental PRL3 inhibitors. Thienopyridone, a potent experimental inhibitor, did not show direct binding to PRL3 but instead inhibited phosphatase activity via redox mechanism. NMR also revealed that other experimental inhibitors did not engage PRL3. Thus, there remains a need to identify potent PRL3-directed inhibitors. Meanwhile, molecular modeling revealed a putative druggable site that has not been thoroughly explored before. The current study provides some scaffolds such as candesartan and particularly, olsalazine, the only binder identified, that could be the starting point of further drug discovery efforts, as well as a putative site that can be targeted in silico. MDMX, a negative regulator of p53, is another important therapeutic target in cancer, along with the homologous protein, MDM2. Inhibitors that block the MDM2-p53 interaction have been identified and despite similarities in the binding site of these homologous proteins, these inhibitors are ineffective against MDMX. It is hypothesized that the flexibility of MDMX contributes to this significant difference in response to inhibitors, despite comparable affinity to their endogenous target, p53. Examination of available inhibitor-bound structures of MDMX reveal a conserved pharmacophore but the structures adopt distinct conformations away from the binding site. This implies that global motions of the protein might contribute to molecular recognition. The conformational heterogeneity in MDMX was further confirmed by collecting residual dipolar couplings (RDCs). Further investigations on both MDMX and MDM2 are necessary to uncover whether the flexibility of MDMX contributes to the differential binding to inhibitors. Finally, NMR relaxation methods and state-of-the-art high-power Carr-Purcell-Meiboom Gill (CPMG) relaxation dispersion measurements, the first documented application on an enzyme, were used to characterize the solution dynamics of GMPK and the changes in dynamics upon GMP binding. Substrate binding resulted in restricting the amplitudes of motion for backbone amide bonds within the picosecond-nanosecond timescale. Meanwhile, CPMG showed dispersion in both in the absence and presence of GMP, such that substrate binding did not quench dynamics within the microsecond-millisecond timescale. Interestingly, more residues are observed to have dispersion in the bound form, some near the C-terminal of helix 3, which has previously been proposed to be involved in product release. Current studies show that substrate binding affect different timescales of protein motion. Future work shall follow how motions within different timescales are affected as GMPK processes its substrates – such as, for instance, binding of ATP analogs within the ATP binding site or simultaneous occupancy of both substrate binding pockets. This paves the way for a complete picture of the relationship of function and dynamics in the conformational enzymatic cycle of a bi-substrate enzyme using GMPK as a model. The current work illustrates some of the diverse applications of NMR on three unique systems that are also drug targets. Information collected here can be leveraged on future structure and dynamics studies as well as drug discovery efforts targeting any of these proteins

A Screen of FDA-Approved Drugs Identifies Inhibitors of Protein Tyrosine Phosphatase 4A3 (PTP4A3 or PRL-3)

Protein tyrosine phosphatase 4A3 (PTP4A3 or PRL-3) is highly expressed in a variety of cancers, where it promotes tumor cell migration and metastasis leading to poor prognosis. Despite its clinical significance, small molecule inhibitors of PRL-3 are lacking. Here, we screened 1443 FDA-approved drugs for their ability to inhibit the activity of the PRL phosphatase family. We identified five specific inhibitors for PRL-3 as well as one selective inhibitor of PRL-2. Additionally, we found nine drugs that broadly and significantly suppressed PRL activity. Two of these broad-spectrum PRL inhibitors, Salirasib and Candesartan, blocked PRL-3-induced migration in human embryonic kidney cells with no impact on cell viability. Both drugs prevented migration of human colorectal cancer cells in a PRL-3 dependent manner and were selective towards PRLs over other phosphatases. In silico modeling revealed that Salirasib binds a putative allosteric site near the WPD loop of PRL-3, while Candesartan binds a potentially novel targetable site adjacent to the CX5R motif. Inhibitor binding at either of these sites is predicted to trap PRL-3 in a closed conformation, preventing substrate binding and inhibiting function

Molecular Dynamics Simulation of an ι-Carrageenan Hexamer as Single and Double Helices

Carrageenan gelation is characterized by a coil-to-helix transition. The conformation the polymer adopts in solution, however, is not yet fully understood at the molecular level. Here, molecular dynamics (MD) simulations were carried out on two possible conformations of ι-carrageenan hexamers – a single helix and a double helix. All-atom simulations were carried out using the Carbohydrate Solution Force Field (CSFF) in explicit water (SPC). Dihedral angle distributions from simulations of the single helix and the individual chains of the double helix were compared revealing similarity between these conformations. In the case of the double helix, an overall attractive interaction has been calculated between the two chains, possibly due to contribution from sulfate groups and hydroxyl groups, preventing it from unravelling. Simulation calculations revealed no significant difference in the total energy of the single helical and double helical carrageenan in equivalent systems indicating no strong preference by the polymer to adopt one configuration over the other. It is suggested that the double helix can form in solution due to the dihedral angle conformation within the single helices and is stabilized by attractive interactions identified. An increased understanding of gelation process, including the coil-to-helix transition, in carrageenans can allow for the development of processing methods to control their properties for specific applications